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NOMENCLATURE 


real  part  of  m 

imaginary  part  of  m 

parameter  in  Eq.  (22) 

derivative  is  the  direction  of  the 
outward  normal  to  surface 

y = 0 (x)  means  lim  y/x  is  bounded 
or  zero  * _ o 


constant  m j 

half-length  of  body  n»2 

horizontal  displacement  due  to  vibration  n 

of  a point  on  body  surface  d/da 

radius  of  body  at  it  • 0 

function  defining  meridian  profile  O(x) 

a constant  in  equation  for  streamline 
body  family  (Eq.  (55)  ) Ps(*)  Legendre  polynomial  of  order  s 

bj/b  Ps  function  defined  by  Eq.  (50) 

non-dimensional  pressure  coefficients  p farfield  pressure 

(Eq.  (58)  ) 

Pq  pressure  in  undisturbed  fluid 

defined  as  C , C for  cos  0 *=  1 n . ........ 

c * 3 Qs  function  defined  by  Eq.  (51) 

the  speed  of- sound  for  the  fluid  _ 7 — 7 : 

r R s W > )J  j.  f/-  t. 


maximum  diameter  of  body 
force  normal  to  surface  of  body 
i,  « 1,  2,  3 components  of  force 
r — b(x*  ) see  Eq.  (23) 


surface  of  body 
u b2  (xp 

subscript  giving  order  of  Legendre 
polynomial  or  Bessel  function 

time 

velocity  of  uniform  streams 
ka  cos  0. 


defined  by  Eq.  (65) 
defined  by  Eq.  (66) 


spherical  Bessel  functions  of  order  s 

2n/\ 

1//32  k 
length  of  body 

i *=  1,  2,  3 direction  cosines 


transverse  velocity  of  a section  of  body 

constant  in  Eq.  (60) 

constant  in  Eq.  (60) 

fluid  velocity  normal  to  surface 

components  of  velocity  (Eq.  (25)  ) 

components  of  velocity 


source  strength 

-U0S'(x,) 


NOMENCLATURE  (Cont'd) 


Xj.Xj.ij  rectangular  coordinates 
*1  = xi  ~ iao<xi>  e'ikct  (Eq.  (23)  ) 

*1  ‘ V* 

/S2  - 1 -M2 

l average  farfield  intensity 


e 


A 

Pn 


M*,> 


P 0 


error  (Eq.  (15)  ) 

cylindrical  coordinate  (Eq  (23)  and 
Fig.  1) 

cos-i  *,/R0 

cos‘1  Xj/Rq 

cos-1  *3/R0 

wavelength 

strength  of  a doublet  with  axis  in 
a given  direction  n 

i *»  1,  2,  3 strength  per  unit  length  of 
the  ith  component  of  the  doublet 
distribution  (Eq.  (41)  ) 

i = 1,  2,  3 coordinates  of  source 
or  doublet 

mass  density  in  undisturbed  fluid 


<P  velocity  potential  of  the  disturbance 

caused  by  the  body 

*T  =-U0*i  + 4. 

*t>s  defined  by  Eq.  (3) 

4>D  defined  by  Eq.  (4) 

4>0  defined  by  Eq.  (24) 

is  such  that  <j>  * <f>  e'ikct  (Eq.  (5)  ) 

Q « cos-  * x j/r 

flj,  flj  constants  in  Eq.  (60) 

( )'  = d(  )/dx, 
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LOW  FREQUENCY  SOUND  RADIATION  FROM  SLENDER 
BODIES  OF  REVOLUTION 


INTRODUCTION 

Since  the  low  frequency  sound  produced  by  surface  ships  and  sub- 
marines can  propagate  over  long  distances  with  small  loss  of  energy, 
it  is  of  interest  to  obtain  some  simple  examples  of  the  radiation  field 
of  such  bodies.  In  this  report  the  vibrating  body  is  considered  to  be  in 
an  unbounded  fluid,  and  thus  the  effect  of  the  free  surface  and  bottom 
(or  other  external  boundary)  is  not  considered.  If  the  body  possesses 
a steady  velocity  along  its  longitudinal  axis,  the  radiated  field  for 
harmonic  vibrations  is  given  for  the  linearized  problem  by  a formula 
which  is  similar  to  the  Helmholtz  formula  and  reduces  to  that  formula 
when  the  steady  velocity  is  zero.  The  results  are  restricted  to  the 
farfield  that  is  due  to  the  vibration  of  a body  of  revolution.  The  body 
length  and  the  wavelength  of  the  vibration  are  much  greater  than  the 
maximum  body  diameter,  and  the  steady  speed  of  advance  of  the  body 
is  much  less  than  the  speed  of  sound  in  the  fluid.  Under  these  con- 
ditions the  formula  for  the  radiation  field  reduces  to  an  integral  over 
a distribution  of  acoustic  sources  and  doublets  along  the  axis  of  the 
body. 

For  the  axially  symmetric  longitudinal  vibrations  of  a body  of 
revolution  only  the  source  term  is  important.  For  transverse  vibra- 
tions the  farfield  is  determined  by  a doublet  distribution.  It  is  shown 
that,  for  the  low  frequency  transverse  vibration  of  a slender  body  of 
revolution,  the  force  which  determines  the  doublet  strength  is  con- 
veniently given  by  the  extended  Lagally  theorem  of  incompressible 


hydrodynamics.  This  theorem  determines  the  forces,  moments,  and 
added  masses  for  a body  which  can  be  represented  by  a closed  stream 
surface  surrounding  a distribution  of  singularities  in  a potential  flow. 
The  theorem  was  given  in  1922  by  Lagally  for  steady  flow.1  Then, 
in  1953,  it  was  extended  to  non-steady  flow  by  Cummins.2 

Related  discussions  of  the  low  frequency  sound  radiation  from 
slender  bodies  have  been  presented  by  Strasberg3  and  Chertock.4*3 
However,  their  examples  are  restricted  to  the  spheroid.  For  the 
spheroid  the  results  given  by  Strasberg  agree  with  those  of  the 
present  report.  Strasberg,  in  turn,  showed  that  his  results  were  in 
agreement  with  those  of  Chertock  for  slender  spheroids. 

If  a given  streamline  body  of  revolution  does  not  differ  too  much 
from  a spheroid,  the  known  results  for  a spheroid  could  be  used  as 
a first  estimate  of  the  radiated  field  of  the  streamline  body.  However, 
the  fact  that  the  streamline  body  is  not  symmetric  about  its  midship 
section  introduces  effects  which  are  not  predicted  from  a spheroid. 

To  illustrate  this,  examples  are  given  of  the  radiated  fields  of  a 
simple  class  of  streamline  bodies  and  of  spheroids  with  the  same 
length  and  diameter. 


1 M.  Lagally,  "Berechnung  der  Krafte  und  Momente  die  stromende  Flussigkeiten  auf 
ihie  Begrenzung  auiiiben,"  Z.  angew.  Math.  Mech.  2,  (1922)  pp.  409-422. 

2W.  E.  Cummins,  "The  Forces  and  Moment  on  a Body  in  a Time-Varying  Potential 
Flow, " Journal  of  Ship  Research  1,  No.  1,  7-18  (1957).  (First  published  as  David  Taylor 
Model  Basin  Report  No.  780,  June  1953.) 

3M.  Strasberg,  "Sound  Radiation  from  Slender  Bodies  in  Axisymmetric  Vibration, " 
Paper  0-28  in  Proceedings  of  the  4th  International  Congress  on  Acoustics,  1962,  Copenhagen 
(Organization  Committee  of  the  4th  ICA  and  Harlang  & Toksvig,  Copenhagen,  1962). 

4G.  Chertock,  "Effects  of  Underwater  Explosions  on  Elastic  Structures,"  Fourth 
Symposium  on  Naval  Hydrodynamics,  B.  L.  Silverstein,  Ed.  (Office  of  Naval  Research 
ACR-92,  1962)  pp.  933-945. 

*G.  Chertock,  "Sound  Radiation  from  Prolate  Spheroids, " J.  Acoust.  Soc.  Am  33 
(1961)  pp.  871-876.  ‘ ’ 
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I.  THE  HELMHOLTZ  FORMULA 


! W 


Consider  the  flow  of  a uniform  stream  about  a slender  body  of 
revolution.  In  Fig.  1,  the  uniform  - U0  is  in  the  negative 


{-a,  o,oj 


Fig.  1 - Coordinate  System 

Xj  -direction,  and  the  body  is  placed  with  its  longitudinal  axis  on  the 
x-axis.  For  a compressible  fluid,  the  velocity  potential  for  the  flow 
about  such  a body  can  be  written  as 

<t>T  = -U0x1  + <K,  (1) 


where 


Ufl  Xj  is  the  velocity  potential  of  the  uniform  stream  and 

$ is  the  velocity  potential  of  the  disturbance  caused  by  the  body. 

(The  disturbance  velocities  are  Wj , i *=  1,  2,  3 with  Wj  = + dQ/dx^ . ) 


Then,  for  small  vibrations  of  the  body,  the  linearized  equation  for 

<X>  is® 


d2<t>  d2Q  «?2  «D 

+ + + 2 — r — — 

dx\  dx2  c2 


d2$> 


C2  <9t2 


(2 


where 

c is  the  speed  of  sound  and 
t is  time. 

It  will  be  assumed  that  the  effect  of  the  vibrating  body  can  be 
approximated  by  changes  in  the  flow  caused  by  distributions  of 
time-varying  sources  and  doublets.  The  solutions  of  Eq.  (2)  for 
a source  and  a doublet  are,  respectively,  the  real  parts  of  7,8 


® 2 * ^3  ^ 


■ ik c ( t - D/c) 


4trR 


(3 


6 1.  E.  Garrick,  "Nonsteady  Wing  Characteristics,"  Section  F of  Aerodynamic 
Components  of  Aircraft  at  High  Speeds,  A.  F.  Anderson  and  H.  R.  Lawrence,  Eds, 
(Princeton  University  Press,  Princeton)  p.  662. 

Reference  6,  pp.  674-675. 

*1.  E.  Garrick,  "On  Moving  Sources  in  Nonsteady  Aerodynamics  and  in  Kirchhoff's 
Formula,"  Proceedings  of  the  First  U.  S.  National  Congress  of  Applied  Mechanics  (The 
American  Society  of  Mechanical  Engineers,  1952)  p.  735. 


is  the  velocity  potential  at  point  (xJt  x2,  x})  and  time  t due  to 
a source  (doublet)  at  £, , £2  > f 3 ; 

is  the  source  strength  (a  source  is  a point  of  outward  radial  flow; 
its  strength  is  the  volume  emitted  per  unit  time); 


is  the  strength  of  .he  ddlblet  with  axis  in  a given  direction  n 


are  the  direction  cosines  of  the  axis  of  the  doublet; 


is  the  wavelength;  and 


9H.  Lamb,  The  Dynamical  Theory  of  Sound  (Edward  Arnold  Lid. , London,  1925, 
2nd  ed. , Dover  Publications,  New  York,  1960)  p.  230. 

10ln  general,  m and  n are  complex  in  Eqs.  (3)  and  (4). 


For  harmonic  time  variation,  let 


<K(xj,  x2.  *3. = 0(xi > *2 * *3 ) c 


Then,  a representation  for  <t>  in  the  region  external  to  the  body,  in 
terms  of  <£  and  its  normal  derivative  on  the  body  surface,  is  given 
by  the  following  generalization  of  the  Helmholtz  formula11 


*2>  x3> 


ii 


where 


d/cln  is  the  derivative  in  the  direction  of  the  outward  normal  to  S. 


The  force  normal  to  ds  is 


Fne  *kCt  ds  = (P  - Pq)  ds  > 


and  the  pressure  p is  given  by 


P -Po  = +*Pokc  <t  + Pouo  T — ’ 

dlj 


11  Reference  8,  p.  737. 


Li*  if.  to  ' 


is  the  pressure  in  the  undisturbed  stream  and 
is  the  mass  density  of  the  fluid. 


Thus 


can  be  written  in  the  form 


For  axialJy  symmetric  flow,  the  surface  integrals  in  Eq.  (12) 
can  be  approximated  by  single  integrals  along  the  axis  of  the  body. 
The  volume  of  fluid  flow  per  unit  time  across  an  elementary  ring 
area  of  the  surface  at  the  mean  position  of  the  body  can  be  written 


The  resultant  of  the  force  F„ 


axis  and  can  be  written  as 


For  points  far  from  the  body  (i.  e.  , terms  proportional  to  1/R' 
are  neglected),  the  following  approximations  hold: 


With  the  conditions  that  M<<1  and  K I M <<1  the  velocity 
potential  <I>  can  be  approximated  by 


where 


With  this  expression  for  <J> , the  pressure  in  the  farfield  is  approxi 
mated  (for  M<  < 1 ) by  the  real  part  of 


A similar  procedure  can  be  applied  for  transverse  vibrations. 
However,  in  Section  Ilia  different  (but  equivalent)  method  is  used 


11.  LONGITUDINAL  VIBRATION 


The  farfield  pressure  is  determined  by  Eqs.  (17)  and  (18)  in 
terms  of  the  flow  at  the  surface  of  the  body.  Let  the  equation  for  the 
undisturbed  body  of  revolution  be 

r = b(a, ) . (19) 

For  the  axially  symmetric  longitudinal  vibration,  the  horizontal 
displacement  of  a point  on  the  body  surface  is  given  by  the  real 
part  of 


ao(x|.') 


1 a, 


0<X1> 


-ikct 


(20) 


The  special  cases  considered  are: 

(a)  rigid  bouy  motion,  aQ  (ij ) = -a0,a0<<b, 

(b)  accordion  motion,  a0(x1 ) = a0  sin 

The  equation  for  the  vibrating  body  is 

f(xr  0,  r,  t)  = r - b(x*  ) = 0, 


(21) 


(22) 


(23) 


where 

x*  = x,  - i a0(x, ) e'lkct  and 


x(,  0,  r are  the  cylindrical  coordinates  shown  in  Fig.  1. 


mm 


* 


For  the  linearized  problem  with  harmonic  time  variation,  the 
axially  symmetric  velocity  potential  will  be  written  as 


♦(*!•  t,  t) 


4>0(x, , r>  + 0(*i> 


r)  e 


-ikct 


(24) 


The  boundary  condition  at  the  body  surface  is 


21 

Dt 


df  1 d( 

w — + w — 
r dr  » r ae 


(25) 


where 

wf  = + d*!*/ dr  , 

Wg  *=  + (?*!>/ <?0  * 0,  for  axially  symmetric  flow,  and 
w j = -UQ  + <?<t>/«9x  j . 

From  Eqs.  (19)  through  (25),  the  boundary  conditions  for  <J>0 
and  <t>  at  the  mean  position  of  the  body  are  (neglecting  second 
order  terms) 


= -Uob'U,), 


(26) 


d 

<9r 


[+iU„  a0'(Xj) 


k c afl  (x  , )J  b'(x, ) , 


(27) 


1 1 


■ 


iJL' 


J 


where 


b'U.)  = 


db(x,) 


d W 


The  boundary  conditions  can  be  applied  at  the  mean  position  of  the 
body, provided  the  slope  b'Otj)  is  of  the  order  of  d/l  «i  and  the 
radius  of  curvature  of  the  body  profile  is  large  compared  with  the 
maximum  amplitude  of  vibration. 

For  rigid  body  motion,  %(*,)  = 0,  and  the  linearized  boundary 
condition  given  by  Eq.  (27)  does  not  depend  on  U0  . For  accordion 
motion,  the  relative  importance  of  the  two  terms  on  the  right-hand 
side  of  Eq.  (27)  depends  essentially  on  the  ratio12 


u n 0 

U.n 


1 u X 
— n M — . 

2 I 


Thus,  for  n = 1,  m = 0.  01,  A/*  = 1 , the  effect  of  the  first  term 

on  the  right-hand  side  of  Eq.  (27)  will  be  negligible.  Only  for 
wavelengths  much  greater  than  the  body  length,  combined  possibly 
with  a large  n,  would  this  term  be  significant  when  M < < 1. 

The  velocity  potential  <t>0  that  satisfies  the  boundary  condition 
(26)  is  determined  approximately  for  a slender  body  by  a distribution 

12See  also  Eq.  (53). 


of  sources,  along  the  axis  of  the  body,  with  strength  per  unit  length 
given  by13 


where 


S(ij ) is  the  sectional  area  of  the  body 


It  will  be  assumed  that  the  velocity  potential 
approximated  at  the  body  surface  by 


where 


This  can  be  approximated  by 


,3S.  Goldstein,  lectures  on  Fluid  Mechanics  (Interscienct  Publishers,  Inc.,  New  York 
1960),p.  183. 


where 


kj  « l/ft2  k,  and 

y 

y = 0(x)  means  lim  — is  bounded  or  zero, 
x -0  * 


The  approximation  requires  that 


kj  tM 


(32) 


Thus,  it  is  not  necessary  for  the  wavelength  A to  be  much  greater 
than  the  body  length  l , provided  M = U0/c  is  sufficiently  small. 
The  radial  velocity  is 


d<6(Xj , r) 
dr 


m(fj)  [l  + ikjM(Xj  - £j)] 


cos  k,  Rj 
(kj  Rl>3 


sin  kjRj 

r cos  k j R j 

sin  k j R j "I 

<klRl>2 

L (kj  Rj  )2 

(k,Rj)}  J 

(33) 


If  the  factor  in  braces  is  expanded  in  terms  of  k,  Rj  then  the 
asymptotic  expansion  for  d<j}/dt  as  r tends  to  zero  is14 


dr  [_  ^2  X^J2,  r 


(34) 


14  Reference  13,  p.  186. 
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' 'k.  : 


Then,  from  Eq.  (27) 


m j(x  j ) = - tca0(i| ) S'(*j ) , 


(35) 


and 


m2(*i)  = + U0aQ  (Jtj)S  '(i,)  , 


(36) 


where 

m(Xj)  = Oj(ij)  + in2(X|)  and 
S(x, ) = irb2  (xj ) . 


(37) 


The  asymptotic  expansions  used  to  determine  d<fj/d r for  small  r 
are  not  valid  at  the  ends  of  the  body  unless  the  ends  are  cusped.1’ 
For  the  vibrating  body,  it  will  be  assumed  that  the  volume  changes 
in  small  regions  near  the  ends  and  the  reaction  on  the  fluid  of  these 
regions  have  a negligible  effect  on  the  farfield  pressure  compared 
with  the  effects  due  to  the  rest  of  the  vibrating  body. 

Since  the  source  distribution  m (x, ) specified  by  Eqs.  (35) 
and  (36)  determines  the  flow  across  the  mean  surface  of  the  vibra- 
ting body,  it  is  also  the  source  strength  for  Eq.  (17)  for  the  far- 
field pressure.  With  this  source  strength,  the  velocity  potential 
t)e-ik‘  1 , as  determined  by  Eq.  (30)  , reduces  to  just  the  first 
term  on  the  right-hand  side  of  Eq.  (17)  (under  the  same  conditions 


15  Reference  13,  p.  185. 


used  in  determining  Eq.  (17)  ).  In  fact,  within  the  approximations 
made,  the  source  distribution  completely  determines  the  farfield. 
There  is  no  net  radial  force  since  the  flow  is  axially  symmetric. 
Hence,  the  last  two  terms  on  the  right-hand  side  of  Eq.  (17)  are 
zero.  The  xi  component  of  force  on  a ring  element  of  area 
2tt  b(x1)dx,  of  the  mean  surface  is  per  unit  length 


Fi<*i>e 


■ike  I 


2 w b (x  j ) 


b'(xj), 

t«b(ij) 


(38) 


where 


[ 


P'Pq  = Po  ikc  r)  +U0 


d<f> 


(*!.')*] 

e-ikc‘. 

j 


It  can  be  verified  that  l/pQ  c F,  (xx)  can  be  neglected  in  Eq.  (17)  in 
comparison  with  m(xt). 

III.  TRANSVERSE  VIBRATION 

For  the  transverse  vibration  of  a body  of  revolution  in  a uniform 
stream,  the  solution  of  Eq.  (2)  is  the  sum  of  the  velocity  potential 
due  to  the  uniform  stream  and  that  due  to  the  vibration.  Let  the 
transverse  vibration  be  such  that  the  circular  sections  of  the  body  do 
not  change  shape  or  radius  but  vibrate  with  velocity  + W (x,)  e'ikct 
perpendicular  to  the  longitudinal  axis.  The  flow  near  the  body  will 
be  considered  as  represented  by  the  superpositon  of  a transverse 
flow  of  velocity  - W (x1)e■‘l<c,  on  the  flow  of  the  uniform  stream 

16 


past  the  undisturbed  body.  Since  the  perturbation  must  be  small  if 
the  problem  is  to  be  linear,  it  will  be  assumed  that  the  maximum 
amplitude  6 of  the  vibration  is  such  that  6 < < b and  the  change  in 
shape  of  the  meridian  profile  is  continuous  and  of  order  0 ( 6 /f  ). 
Then,  with  the  additional  conditions  that  d/I  , r/\,  U0/c,  W/c  < < 1, 
the  solution  of  Eq.  (2)  that  holds  near  the  surface  of  the  body  is  16 

*<*,.'•  n)e*ikct  = _“S(Xj)W(*j)  — - e-»kct  [ i + o(k2rJ  In  kr ) ] , (39) 

where 

W (x  j ) >s  the  transverse  velocity  of  a section  of  the  body  and 
X3 

cos  Q ~ — . 

r 

For  incompressible  flow  this  same  solution  holds  with  error  factor 
[l  + 0 (d2/l  2 In  d/I ) ] ,17  Hence,  the  flow  can  be  considered  as  quasi- 
stationary (i.  e. , time-varying-incompressible)  flow  at  the  surface 
of  the  body.  It  can  be  represented  by  a distribution  of  three-dimen- 
sional doublets  with  axes  perpendicular  to  the  body  axis  and  in  the 
direction  of  w (x,)  . The  strength  of  the  doublet  distribution  per 
unit  length  along  the  axis  is 

*■}(*,,  0 = + 2 S(x  j ) W(x  j ) e'* k c * . (40) 


16J.  W.  Mile*,  "On  Non-Steady  Motion  of  Slender  Bodies,  " Aeronaut.  Quart.  2, 
November  (1950)  p.  186. 

17B.  Thwaite*,  Incompressible  Aerodynamics  (Oxford.  Clarendon  Press,  1960)  p.  393. 


For  the  flow  of  the  uniform  stream  past  the  body  there  is  a doublet 
distribution  p.j  (xt)  along  the  axis  with  strength 

/*,(*,)  = U0S<*,).  (41) 


However,  with  the  conditions  imposed,  the  approximate  solution  of 
Eq.  (2),  as  given  by  Eq.  (38)  , does  not  depend  on  U0,  and  the 
reaction  of  the  body  on  the  fluid  is  approximated  by  the  reaction  of 
the  time-varying  doublet  distribution  pj . Applying  the  extended 
Lagally  theorem  to  this  doublet  distribution  gives  the  reaction18’19 

FJ(«1,t)d*i  = ~P0-|;  Mxi>  t)d*i-  (42) 

= +2  i p0  Icc  S(*j ) W(x, ) e',kcl  dxj  . (43) 

The  mutual  interactions  of  the  doublets  would  involve  products  of 
the  velocities  at  different  points  on  the  body  and  are  neglected  in  the 
linear  theory. 

For  the  slender  body  the  force  F3  dx,  exerted  on  the  fluid  at 
each  section  acts  as  an  isolated  doublet  with  axis  in  the  x3 -direction. 
The  resulting  farfield  velocity  potential  is 


<t>(x, , x2,  x},  t)  = -■ 


.-ike  (i  • D<)/c) 

4 TT  R 


cos  #3 \ 2 i 


ik  J8  S(^)W(^)e-ik  cos  6°,,fl  df,|  (44) 


where 


18  Reference  2,  p.  18. 

19 L.  Landweber  and  C.  S.  Yih,  "Forces,  Moments,  and  Added  Masses  for  Rankine 
Bodies,"].  Fluid  Mech.  1,  (1956)  p.  332. 
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The  farfield  pressure  is  given  by 


P - P0  = +T-  P0  k2c  cos  0 3 


-ilcc  (i  - Dq/c)  r a 


f“  W(^)S(^)e-ikcos  d£,. 

•/-a 


The  same  result  is  obtained  if  the  velocity  potential  given  by 
Eq.  (39)  is  used  in  Eq.  (12).  In  this  case  the  two  terms  of  Eq.  (12) 
are  equal.  The  first  term  gives  the  effect  of  the  acceleration  of  the 
displaced  fluid,and  the  second  term  gives  the  added  mass  effect  in 
the  approximating  two-dimensional  incompressible  flow  near  the 
body. 


IV.  EXAMPLES  OF  AXIALLY  SYMMETRIC  VIBRATION 
For  axially  symmetric  vibration,  the  farfield  pressure  is  given 


• -ilcc  (t  -Do/ c)  r a 

P-Po  = ~ 7“  Po  k c 5 I 

’ 77  *'0  J- a 


cos  9q  • £ j 


where  Eqs.  (17)  and  (18)  were  used  and  the  force  term  of  Eq.  (17) 
was  neglected  in  accordance  with  the  analysis  indicated  at  the  end 
of  Section  II.  The  types  of  vibration  are  specified  by  Eqs.  (20)  and 
(21).  Using  Eqs.  (35)  and  (36)  yields  the  corresponding  source 


strengths: 


(a)  rigid  body  motion,  m j (x  , ) = + k c a0  S '(x  t ) , m 2 (x  , ) =0, 


(b)  accordion  motion,  m , (x  , ) ■=  -kca 


0 sin  (t^)  s (Xl)  ’ 
•*<«,>  -7X7-^) 


(48) 


where 

m(xj)  is  defined  by  Eq.  (37), 


and 


S'(x,) 


dS(x, ) 


dr, 


It  will  be  assumed  that  the  slope  of  the  sectional  area  curve  can  be 
represented  by  a finite  sum  of  Legendre  polynomials. 20  Thus, 


S-U,)  -£*.!>. 


s = 0 


£)■ 


(49) 


where 

the  A.  are  constants. 


With  the  source  strength  defined  by  Eqs,  (48)  and  (49),  the  integral 
in  Eq.  (46)  is  the  sum  of  integrals  of  the  form 


= i(-i)s  a I js  (u  + v ) - j (u  - v ) | , 


20  For  the  spheroid  and  streamline  body  examples  considered  in  this  section,  S '(x  ,) 
is  given  exactly  by  (4!))  with  only  a few  terms.  For  mote  general  forms,  S/(x()  can  be 


approximated  by  (49), 


(50) 


A - 
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and 


g>  ‘ £ (t  t)  f.  (t)  80 ' ' f ' “ (. 


= (-i  )s  a l j (u  + v ) + j (u  - v)  J , 


(51) 


where 


u = It  a cos  dQ  , 
v = n rr/2  , and 

j ( u)  is  the  spherical  Bessel  (unction  of  order  s. 
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With  these  results,  the  farfield  pressure  for  rigid  body  motion  and 
accordion  motion  are,  respectively,  the  real  parts  of 


i e-*kc  (t-  D0/c ) _s 

P " Po  = ~ Po  k c a ao  z 

1 V K0  sTO 


£ (-ns  as  js(u) , 


(52) 


and 


1 .22  e 

P - Po  = - — P„  k c a an  - 


ike  (t  - D 0/c  ) s 


R 


0 s =0 


Y ( ~ i)S  A s 1 1 js  (u  + v)  - js  (u  - v] 


(53) 


-"My  ( js  (u  + v)  + js  (u  - v)l  | . 


2IP.  M.  Morse  and  H.  Fcshback,  Methods  of  Theoretical  Physics  (McGraw-Hill  Book 
Co.,  Inc.,  New  York,  1953)^.  1573. 


21 


»»  ** 


•t  t,  Stamms 
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In  the  following  examples  it  will  be  assumed  that  nM(A/l)<  < 1,  and 
the  second  terms  of  Eq.  (53)  will  be  neglected. 

The  equation  of  the  meridian  profile  of  a spheroid  with  center  at 
the  origin  and  axis  of  revolution  along  the  x,  -axis  is 


r(x , ) - b yl  t , 


where 


r2  = *2  + x2. 


a is  the  semi-major  axis, 


b is  the  semi-minor  axis,  and 


x,  = Xj/a  . 


To  illustrate  the  effect  of  a change  in  the  distribution  of  volume 
along  the  axis  of  revolution,  a simple  family  of  streamline  bodies 
can  be  constructed  by  adding  the  following  function  to  the  equation 
for  the  meridian  profile  of  a spheroid. 

r | (x  | ) = 2bb,x,  yjl  - x f , (55) 

where 

b,  = bj/b  and 

bj  is  a constant  that  determines  the  amount  of  distortion  from  a spheroid. 

Figure  2 shows  the  meridian  profile  of  a spheroid  with  b/a  = 0.  1 
and  the  profiles  of  two  streamline  bodies  with  b/a  = 0.  1,  bt  = 0.  3 
and  b/a  = 0.  1,  b,  = 0.  5. 


■B 


22 


The  slope  of  the  sectional  area  curves  for  the  spheroid  and  the 
streamline  bodies  can  be  written  as 


S'(x,) 


-2  jt  — P . <T,  ) 
a 1 1 


(spheroid)  , 


(56) 


and 


S'(x,)  = 


P.fiT,)  +8  b1P2(5'1) 


(streamline  body). 


(57) 


The  real  part  of  the  farfield  pressure  determined  by  Eqs.  (52)  or 
(53)  can  be  written  as 


<P  -Po>  . = P0  k2c2b2 


real 


*0  |CS 


sin  (kc(t  - D0/ c ) I 


+ C, 


cos  [kc(t  - D0/c)  ] 


}■ 


(58) 


where 

C$  and  Cc  are  non-dimensional  functions  of  the  parameters 
specifying  the  body  and  motion. 
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t V. 


Table  1 

E FUNCTIONS  Cs  AND  C 


For  definition  of  J*  and  J*  see  Eqs.  (65)  and  (66) 


Formulas  for  Cs  and  Cc  for  the  spheroid  and  streamline  body 
families  with  rigid  body  and  accordion  motion  are  given  in  Table  1. 
For  the  spheroid  in  rigid  body  vibration  and  the  simplest  case 
(i.  e.  , n = 1)  of  accordion  vibration  these  results  agree  with  those 
ot  Strasberg,J  who  has  shown  that  these  results  agree  with  those 
given  by  Chertock5  for  slender  spheroids  and  u < 3.  Since 
Chertock's  results  are  based  on  an  asymptotic  expansion  (for  b/l  - 0) 
of  the  exact  solution  for  a spheroid,  the  condition  u < 3 represents 
a restriction  on  the  results  based  on  the  slender  body  theory  of  the 
present  report. 

For  a given  body  and  mode  of  vibration,  the  independent  vari- 
able in  the  formulas  for  Cs  and  Cc  is  u = ka  cos  0O.  Thus, 
with  ka  held  constant,  the  formulas  show  how  the  axially  sym- 
metric farfield  pressure  given  by  Eq.  (58)  depends  on  the  position 
of  the  observation  point  as  specified  by  cos  0O  = x,/!^  . For  the 
spheroid  in  accordion  motion  both  the  body  and  the  motion  are 
symmetric  with  respect  to  xp  and  only  a symmetric  formula 
occurs.  Similarly,  for  rigid  body  motion,  which  is  not  symmetric 
with  respect  to  xi,  only  a skewsymmetric  formula  occurs  for  the 
spheroid.  The  meridian  profile  of  the  streamline  body  is  not 
symmetric  with  respect  to  xIt  and  as  a result  one  of  the  functions 
Cs  * Cc  is  symmetric  and  one  is  skewsymmetric  with  respect  to 
x j in  each  case.  These  functions  are  shown  in  Figs.  3 and  4.  A 
further  comparison  of  the  radiated  fields  is  provided  if  the  average 
farfield  intensity,  based  on  the  linear  theory,  is  defined  by 
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if  a a ^ , a „ *•  '' 


is  plotted  in  Figs.  5 and  6 


The  function  C 


sPHeaoiD 


sntAMUNi  body  6,0 


STREAMLINE  BODY  6,0 


SPHEROID 


STREAMLINE  BODY  6, *.3 


STREAMLINE  BODY  6, 


Fig.  5 - Intensity  Rigid  Body  Vibration 


Fig.  6 - Intensity  Accordion  Vibration 
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V.  EXAMPLES  OF  TRANSVERSE  VIBRATION 

For  transverse  vibration  the  farfield  pressure  is  determined  by 
Eq.  (45).  Let  the  transverse  velocity  of  the  sections  be 

*(*!.«)-  (60> 

where 

Wl,W2,Sll,  and  fl2  are  constants. 

If  the  sectional  area  of  a given  body  is  expressed  as  a sum  of 
Legendre  polynomials,  Eqs.  (50)  and  (51)  can  be  used  to  evaluate 
the  integrals  appearing  in  the  expression  for  the  farfield  pressure. 
For  the  spheroid  and  the  streamline  bodies  considered  in  Section  IV, 
the  sectional  areas  are,  respectively, 

S(x j ) = >rb2  (l  -If2)  (hi) 


= irb2  — [PQ  - P2(lf1 ) ] , 


S(xj)  - *b2[(l-l2)  +41,(1, -*»»)+ It*)] 

- 'i>!{  juvV-'VV'  - 


The  farfield  pressures  for  the  spheroid  and  the  streamline  bodie: 
are,  respectively, 


P-P0 


+ pQ  k2  c2  b2  cos 


e-ikc  (t  - Dg/ c ) 

*~o 


and 


(.It), 

~7V‘  °‘l * [(i**sj) 


where 


’Jdr)  ■•(-£) 
(“4‘)  + (~£)  ' 


* - 1,  2 ; 


and 


(■•¥)  (-¥-) 


( 


If  W2  = 0 and  ft,  =1,  the  maximum  amplitude  of  vibration  is  at 
the  half-length  of  the  body,  and  there  are  nodes  at  the  ends.  Using 
the  notation  of  Eq.  (58)  with  a0  = - a W,  /c  yields  the  expressions 
for  C,  and  Cc  given  in  Table  1.  For  W,  =0,  ft2  = 1»  the 
maximum  amplitude  of  vibration  occurs  at  the  ends  of  the  body, 
and  there  is  a node  at  the  half-length.  Setting  a<>  = - a W2  /c,  the 
expressions  for  C,  and  Cc  are  given  in  the  last  line  of  Table  1. 

(V 

Figures  7 and  8 show  the  functions  Cs  and  Cc  (where  C, 

= C,  cos  83  and  Cc  = Cc  cos  83  ) for  these  examples  of  trans- 
verse vibration  for  spheroid  and  streamline  body  families.  The 
corresponding  intensities  are  given  in  Figs.  9 and  10. 
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Fig.  7 - C*  and  Cc  Transverse  Vibration 
w2  ■ °,  fl,  ■ 1 
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Fig.  8 - C»  and  Cc  Transverse  Vibration 
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Fig.  9 - Intensity  Transverse  Vibration  Fig.  10  - Intensity  Transverse  Vibration 

w2  * o,  flj  = i Wj  * o,  n2  ■ i 

As  a final  remark  on  the  comparison  shown  in  Figs.  3-10,  it 
should  be  noted  that  the  streamline  body  family,  as  defined  by 
Eqs.  (54)  and  (55),  was  selected  to  give  relatively  simple  expres- 
sions for  the  farfield  pressure.  All  comparisons  are  for  bodies  with 
the  same  length  I = 2a  and  the  same  radius  b at  the  center.  This 
is  illustrated  by  Fig.  2.  It  might  be  expected  that  the  comparison 
should  be  for  the  same  ratio  of  maximum  diameter  to  length,  or  for 
the  same  volume,  or  some  other  specification  of  the  spheroid  to  use 
in  comparison  with  a streamline  body.  However,  consider  Fig.  3, 
which  shows  Cs  and  Cc  for  rigid  body  vibration.  The  spheroid 
provides  a good  approximation  to  Cc  for  the  steamline  bodies, 
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but  no  spheroid  provides  an  approximation  for  C,.  Thus,  no 
significant  improvement  in  the  approximation  of  the  farfield  of  the 
streamline  bodies  would  be  obtained  by  selecting  a different 
spheroid.  Similar  remarks  hold  for  the  other  examples. 
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